Spin-orbit interaction and spin relaxation in a two-dimensional electron gas 
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Using time-resolved Faraday rotation, the drift-induced spin-orbit field of a two-dimensional elec- 
tron gas in an InGaAs quantum well is measured. Including measurements of the electron mobility, 
the Dresselhaus and Rashba coefficients are determined as a function of temperature between 10 
and 80 K. By comparing the relative size of these terms with a measured in-plane anisotropy of the 
spin-dephasing rate, the D'yakonv-Perel' contribution to spin dephasing is estimated. The mea- 
sured dephasing rate is significantly larger than this, which can only partially be explained by an 
inhomogeneous g factor. 



The possibility to manipulate spins in semiconduc- 
tors is a requirement for future spin-based information 
processing.— Using the spin-orbit (SO) interaction^ii^ is 
a promising way to precisely control spin polarization 
because of its simple principle based on external gate 
electrodes4i^ Manipulation of spins using the SO inter- 
action has been shown in various semiconductor systems, 
such as bulk semiconductors^, two-dimensional electron 
gasesi (2DEGs), and even quantum dots containing only 
one single electron^ On the other hand, the SO interac- 
tion is a source for spin dephasing. In 2DEGs, the SO 
interaction induces a linear /c-dependent splitting.^ This 
splitting gives rise to an effective magnetic field, leading 
to dephasing of the polarized electron spins This effect 
is known as the D'yakonv-Perel' (DP) mechanism, and its 
control through manipulation of the SO interaction has 
been proposedii as an alternative to the ballistic spin 
transistor^. A careful engineering of the SO interaction 
is therefore crucial for using it to manipulate the spin. 

In a 2DEG at intermediate temperatures, it is often as- 
sumed that the spin decay is governed by the DP mecha- 
nism. ^^'^■^ Based on this assumption, information on the 
SO interaction in semiconductor quantum wells (QWs) 
was obtained from measurements of the spin-dephasing 
rate li^ii^ii^ An independent measurement of the relative 
size of the SO interaction in (llO)-grown QWs u sing the 
photogalvanic effect has been described in Ref. [13 and 
compared to the spin decay time. In this paper, we re- 
port on quantitative and independent measurements of 
the SO interaction and the spin-dephasing rate in an In- 
GaAs QW, utilizing time-resolved Faraday rotation. In a 
further development of the method described in Ref. [tI, a 
well-defined current is applied in the 2DEG using Ohmic 
contacts and a mesa structure (in Ref. |7|, the electron 
drift was induced by an ac voltage applied to Schottky 
contacts in an unstructured 2DEG). The drifting elec- 
trons see an effective SO magnetic field, in the following 
referred to as drift SO field. The sizes of its two con- 
tributions, the Rashba'^ and the Dresselhaus^ fields, are 
determined as a function of temperature T from the mea- 
sured influence of the in-plane electron drift velocity on 
the spin precession. Comparing our results with mea- 
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FIG. 1: (a) The InGaAs/GaAs QW is shaped into a cross 
with 150-/im-wide arms contacted by Ohmic contacts. Ad- 
ditional contacts on one arm of the cross allow a four-point 
measurement of the voltage drop Vxx and the determination 
of electron sheet density and mobility, (b) 6 and 4> are the 
angles of the magnetic and electric fields with respect to the 
[110] axis. Btot is the vector sum of the external magnetic 
field and the two SO effective magnetic field contributions. 



sured spin-dephasing rates and their in-plane anisotropy, 
we find that DP is not the only mechanism for spin decay 
in our samples at T between 10 and 80 K. 

The 2DEG we use in this work is located in an 
Ino.iGao.gAs/GaAs QW. Electrons are confined to a 20- 
nm-thick Ino.iGao.gAs layer that is n-doped (3x10^^ 
cm^'') to ensure a small electron scattering time such 
that we are in the dirty limit of the SO interaction-, 
where the frequency loso of spin precession about the 
SO fields is small compared with the momentum scat- 
tering rate l/rp (i^so''> ~ 10~^ for our samples). On 
both sides of this layer, there is a 10-nm-thick GaAs 
spacer layer and a 10-nm-thick layer of n-doped GaAs. 
A 10 nm cap of undoped GaAs completes the structure, 
grown by molecular beam epitaxy and forming a 2DEG 
40 nm below the surface. We use wet etching to pattern 
a cross-shaped mesa as shown in Fig. [Tfa) and create 
standard AuGe Ohmic contacts in the four ends of the 
cross. Four additional contacts on one arm of the cross 
allow its use as a Hall bar to determine the resistivity 
and carrier density of the 2DEG. Two samples with the 
same structure are glued into one chip carrier, whereby 
one sample is rotated by 90° to allow the SO interaction 
to be measured in one cool-down. At 40 K, the two-point 
resistance of the crosses in x or y direction is 4.1 kil. 
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We use additional resistors Rs = 4.7 kfi to compensate 
for small variations in the contact resistance and apply 
voltages Vi — Va cos(0) and V2 = Va sm{(j)) as shown in 
Fig. (TJa). All angles are given with respect to the x axis 
along [110], as seen in Fig.llfb). We obtain the resistivity 
of the 2DEG during optical experiments monitoring the 
ratio of the voltage drop V^x and the current Iq through 
one arm of the mesa [see Fig. [IJa)], and measure a value 
of 770 il/sq at 40 K. The voltages Vi and V2 create an 
electric field E in the center of the cross in direction (j) 
and with an amplitude proportional to Va- Because Vxx 
is monitored, the corresponding component of E can be 
determined directly (see below). The electric field shifts 
the Fermi circle by an amount of Sk — m*^E/h, where 
m* is the effective electron mass, /i is the electron mo- 
bility, and h is Planck's constant divided by 2-k. In the 
dirty SO limit, the shift induces drift SO fields that can 
be expressed as^ 
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with 6k ~ {Skx, Sky), g is the electron g factor, iib is the 
Bohr magneton, and a and P are the Rashba and Dres- 
selhaus spin-orbit coefficients, respectively. The Rashba 
field has its origin in the structure inversion asymme- 
try (SIA) due to nonuniform doping on both sides of the 
QW and the presence of the surface on one side of the 
QW. The Dresselhaus field is a consequence of the bulk 
inversion asymmetry (BIA) of the zinc-blende structure. 
Cubic Dresselhaus terms do not change the linearity of 
Bbia in 6k, but introduce a correction of 1 — \kp/{kl), 
where fcp is the Fermi wave number and (fc^) the expecta- 
tion value of the squared wave number along the growth 
direction z. Taking a sheet density of 5.2 x 10^^ m^'^ 
(see below) and approximating (k^) by {n/w)'^, where 
w = 20nm is the QW width, we obtain kp « 1.32(fc^). 
This gives a correction in /3 of about than 35%, which 
will be neglected in the following. 

An external magnetic field i?cxt is applied in the direc- 
tion 6 as seen in Fig. [TJb). The angles 6 are 90° for sam- 
ple 1 and 180° for sample 2. A transverse electron polar- 
ization precesses coherently about the vector su m^^i^^i^° 
of Bcxt and the drift SO fields defined in Eq. ([T|) with a 
frequency given by the modulus of this total field vector. 
If Bext ^ Bsia,Bbia, the total field can be approxi- 
mated asi 
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Because of the difi'erent angular dependencies of the 
Rashba and Dresselhaus SO fields, the two contribu- 
tions can be distinguished. We use time-resolved Fara- 
day rotation to determine the Larmor frequency Ol = 
gUsBtot/f^ of the spins precessing about Btot- For this, 
the output of a pulsed Ti:sapphire laser with a repeti- 
tion rate of 80 MHz is split into a pump and a probe 
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FIG. 2: (a) [Btotj in the center of the cross (symbols) as a 
function of the measured electric field. The slope of the fitted 
line is 2.71x10"*^ Tm/V. (b) Dependence of the fitted slope 
on the position x (crosses). Assuming a constant a — /3 of 
4.9 X 10"^'' eVm, a simulation of the electric field reproduces 
the measured values (solid line). 



beam. The pump/probe intensity ratio is 10/1 and, un- 
less stated otherwise, the pump power is 500 fiW, focused 
onto a spot 30 /im in diameter. The circularly polar- 
ized pump pulse is tuned to the absorption edge of the 
QW at 1.44 eV and creates a spin polarization along the 
growth axis of the QW. With a pump power of 500 /iW 
and assuming an absorption of about l%2i we obtain 
a photo-excited carrier concentration on the order of a 
few lO^'^cm"^, which is more than a magnitude smaller 
than the equilibrium carrier sheet density in the QW (see 
below). The Faraday rotation of the linear polarization 
axis of the probe pulse transmitted is proportional to the 
spin polarization along the QW growth axis. Changing 
the delay At between pump and probe reveals the spin 
dynamics of the system, and the Faraday rotation angle 
can be described by ylexp(— At/r2*) cos(f2LAi). Here, A 
is the amplitude of the Faraday signal and T2 the spin 
dephasing time. A measurement of in a known mag- 
netic field reveals an electron g factor of —0.29, assuming 
that the g factor is negative. 

Figure HJa) shows i?tot measured in the center of the 
cross as a function of the electric field Eq — Vxx /I be- 
tween the two contacts in the right arm. These con- 
tacts are separated by a distance of I = 100 ^m. V2 is 
set to ground. The temperature of the sample is 40 K, 
and Bext is oriented along the [110] direction, therefore 
6 = (j) = 90°. The data in Fig. [D^a) contain values from 
sweeps of Vi up and down. The up and down sweeps fit 
nicely to a straight line, showing that we can exclude a 
drift of over time, which might be caused by nuclear 
polarization or drift of laser power or temperature. 

As a result of the sample geometry and the grounded 
contacts to the left and the right of the center, the electric 
field in the center of the cross is reduced as compared to 
the value Eq measured in the arm. By scanning the laser 
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FIG. 3: Btot as a function of the direction (j) of the in-plane 
electric field for difi'erent amplitudes E and for two external 
magnetic field directions (a) 9 = 180° and (b) 6* = 90°. The 
solid lines are fits using Eq. ([2]| . The data are obtained at 
T = 20K. 



beam along the x-axis and centered in the j/-direction 
[see Fig. HJa)], we obtain a cross section of the drift SO 
field that is related to the electric field distribution.^^ 
The resulting slopes of the linear fits of Stot vs. Eq are 
shown as a function of x in Fig. [H^b). We see a pro- 
nounced dip in the center of the cross, which is explained 
by the reduced electric field. The solid line in Fig. ^h) 
represents the solution of a numeric simulation of the 
electrostatics using a two-dimensional partial differential 
equation solver. Assuming that a and (3 are independent 
of the position and using the measured mobility (see be- 
low), the only fit parameter left is the difference a — /3. 
The measurement and the simulation are consistent and 
show that the electric field E in the middle of the cross 
is 0.71 times smaller than the measured value Eq. This 
correction is taken into account in the following when 
indicating electric-field values. 

To disentangle a and f3, we position the laser spot in 
the center of the cross and rotate the electric field. We did 
such experiments for different amplitudes of the electric 
field up to 3.3kV/m and for two different configurations 
of Bext- In Fig.ISta), e = 180°, and in Fig.[3i;b), 9 = 90°. 
The data are obtained at 20 K. Btot oscillates in (f> with an 
amplitude that is proportional to a + (3 ioi 9 = 180° and 
to a — (3 for 9 = 90°. The difference in the amplitude 
for the two cases (note the different scales) shows that 
-BsiA and -Bbia are comparable in relative strength and 
that the interplay of the two SO effects gives rise to an 
anisotropic spin splitting in k space. The solid lines are 
a fit to the data using Eq. Small deviations of the 
data from theory in the (f) direction could result from a 
slight accidental off-center position of the laser spot. 

To calculate the SO coefficients a and f3 from the mea- 
sured -BsiA and -Bbia, we need a value for the drift mo- 
mentum of the electrons. This is obtained from a Hall 
measurement of the sheet resistivity p and the sheet car- 



rier density n. We calculate the mobility fi using /i = 
1/nep, e being the electron charge. In the dark, the re- 
sistivity of the 2DEG is approximately 1000 O/sq and de- 
creases to 770 ri/sq under conditions of the optical mea- 
surements. The Hall sheet densities are 5.8 x 10^^ m^^ 
under illumination and 4.5 x 10^^ m^^ in the dark. These 
densities are constant from 10 to 80 K. From Shubnikov- 
de Haas oscillations at T < 20 K, we get a sheet density 
of 5.2 X 10^^ m^^ under optical illumination. A small 
parallel conductivity of a doping layer could explain the 
difference in the two numbers. Such a parallel conduc- 
tivity does not influence the optical experiments as these 
electrons do not contribute to the Faraday signal. 

The mobilities extracted from the Hall measurement 
under illumination are shown in the upper inset of 
Fig. Ufa). The mobility does not change significantly over 
the temperature range measured. Using the results of the 
transport measurements and assuming that m* = 0.058, 
we can calculate Sk and use Eq. ([1]) to obtain a and (3 for 
all the temperatures measured. The results are displayed 
in Fig. Hl^a). Error bars show a 2a confidence interval. 
The wafer used for this work is the same as the one used 
for sample 3 in Rcf. 7|. We measure values for a and /? 
that are by a factor of 2-3 smaller, which we attribute to 
a more precise determination of the electric field in this 
work. In addition, different wafer processing and oxida- 
tion of the wafer surface over time might influence the 
SO coefficients measured in this shallow 2DEG. Varia- 
tions in a and (3 for subsequent cool downs are within 
the error bar. Note that we extrapolate the mobility 
and the electric field in the center of the cross from a 
transport measurement away from the center. We can- 
not exclude that we underestimate the absolute values for 
a and P because of a reduced electron drift momentum 
that might result from, e.g., screening by the optically 
excited charge carriers. The model of Viha^'^ with the 
results and the parameters used by Hiibner et ali^ pre- 
dicts the T dependence of the band parameters, from 
which we estimate the T dependence of a and /? using 
k ■ p theoryj^ The calculated T-induced change in a and 
P between 10 and 80 K is in the sub-percentage range 
and thus much smaller than our measurement error. In- 
terestingly, in Ref. a linear increase in a with T for 
higher T was observed on a [110] QW, which the authors 
could not explain with k ■ p theory. 

We find no dependence of the SO coefficients on Bcxt, 
in agreement with our assumption that the precession 
frequency is given by the modulus of the vector sum of 
Bcxt and the drift SO fields given in Eq. (P). Figure Hljb) 
shows a measurement of a — /3 vs i?cxt ■ The insensitivity 
of the result on -Bcxt excludes a significant admixture of 
a k dependent and anisotropic g factor, as was stipulated 
in Ref. 26. To test the reliability of our method, we also 
checked whether a lower pump power will influence the 
outcome of the measurement. This could occur from, 
e.g., a population of higher energy states with larger 
pump power. We found, however, that a — (3 does not 
depend significantly on the pump power, as seen in 
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FIG. 4: (a) Measured fi, a and /? vs. T. (b) Measured a- P 
shows no significant dependence on Bext- (c) Pump-power 
dependence of a — /3. Tfie data in (a) and (b) were obtained 
witli a pump power of 500 /iW, tlie data in (b) and (c) at 
T = 40 K. 



Fig. He). 

When a and f3 are of similar magnitude, the spin hfe- 
time is strongly anisotropic with respect to the direction 
of Boxt in the plane of the 2DEG^iMii6^^^ This 
anisotropy is a consequence of the DP mechanism be- 
cause the spins precess about a SO field whose direction 
becomes independent of k for a ~ (3. From the mea- 
sured anisotropy in Tj* and the relative size of a and /3, 
we estimate the contribution from the DP mechanism to 
the spin dephasing. In Fig. EJa) , the spin relaxation rate 
I/T2* is plotted as a function of T for the two orienta- 
tions of i3oxt=0.99T. The T dependence is rather small, 
and we can clearly see an anisotropy of I/T2*, confirm- 
ing the anisotropic spin splitting in our system. Because 
our 2DEG is well in the dirty SO limit, we can use the 
motional narrowing limit of the DP mechanism^ where 
the spin dephasing due to |fc|-dependent SO fields is de- 
creased by spin-preserving scattering. This gives the fol- 
lowing expressions for the anisotropic spin decay rates: 
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x^y^z are the relaxation times of spins oriented 
along X, y, or 2;||[001]. C is a constant that depends 
on T, Fermi energy, scattering time, and the scattering 
mechanismjSS, If we apply a large external magnetic field 
i^LTx,y,z ^ 1 ), we can write the DP spin-dephasing rate 
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For the difference, we get t^^qj — t-{iso) ~ ^'^(^ read 
a value of about 0.4 x 10^ s~^ in Fig. ^du) for the dif- 
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FIG. 5: (a) Spin depfiasing rate I/T2 for two different in- 
plane directions 6 vs. T from 10 to 80 K. The average of 
three measurements with 50 /iW pump power was used, (b) 
I/T2 vs. Bcxt and linear fits for three different pump powers. 
The data is measured in the 9 = 90° configuration at 10 K, 
and at 40 K for the 500 /iW case, (c) Pump-power dependence 
of 1 /T2 measured at 10 K in the 6 — 90° configuration. 



ference. Using the measured values for a and f3, we get 
C = 6.6 X 1035m~2gv-2s-i. From Eqs. © and Q, this 
yields relaxation rates for DP of about 0.6 x 10^ for 
9 = 180°, and of 1.0 x 10^ s^^ for 9 = 90°. In Eq. ©, 
it is assumed that the SO splitting is linear in k. As 
mentioned earlier, we are in a regime where kp « (^z)- 
Taking into account the cubic Dresselhaus terma^S we 
find only a small correction to the values for the spin 
relaxation rate obtained above. As the total I/T2 lies 
between 2.1x10^ and 2.5x10^ s~^, other spin-dephasing 
mechanism must be present in our sample. 

To exclude optical recombination as a source of de- 
cay of the Faraday signal, we measured the time-resolved 
reflectionji^ which exponentially decays with a decay 
time of less than 100 ps (data not shown). Interpreting 
this time as the electron-hole recombination time pro- 
vides evidence that the spin polarization, which is ob- 
servable over a much longer time scale in the Faraday 
signal, must get imprinted onto the equilibrium electrons 
in the QW conductance band through recombination of 
unpolarized electrons and holes^ It is therefore justified 
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Relaxation mechanism DP Ag Sum Measured 
61 = 90° (l/s)xlO^ 0.6 0.7 1.3±0.2 2.1± 0.1 
6* =180° (l/s)xl0^ 1.0 0.7 1.7±0.2 2.5± 0.1 



TABLE I: Contributions to the spin decay rate I/T2 at 30 K, 
as estimated from the measured anisotropy of I/T2 , SO con- 
stants Q and (3, the _B-dependence of I/T2 , and the electron 
mobility. Measured values are obtained at a pump power of 
50AtW. 

to interpret the decay time of the Faraday signal fitted 
in a window from 80 to 1000 ps as the decay time T2*. 

From the dependence of Tj* ^cxt, information on 
the mechanism of spin dephasing can be obtained. The 
DP spin dephasing rate does not depend on B^xt in the 
motional narrowing regime and for flLTp <C l."^^ In con- 
trast to this, a B-dependence that is intrinsic to the DP 
mechanism is observed in high-mobility samples. In our 
low-mobility samples, we find a linear increase in with 
Bcxt, as shown in Fig. [SJ^b). Such a linear i3-dependence 
is evidence of an inhomogeneous dephasing due to a vari- 
ation Ag of the g factor in the area of the 2DEG probed, 
described by a dephasing rate l/rAgr^^ 

1 _ Ag^sBcxt 



We suspect that the g factor variation could be a conse- 
quence of the in-well doping. That the sample is rather 
inhomogeneous is also seen in a photoluminescence ex- 
periment, in which we observe a broad luminescence peak 
from the QW (not shown) with a full width at half max- 
imum of about 20meV. We used different pump intensi- 
ties and find a similar slope for the dashed linear fits in 
Fig. IHb). The 500 /iW measurement was done at 40 K, 
and the other two at 10 K. From these data, we con- 
clude that Ag is quite constant for different pump pow- 
ers. From the slopes in Fig. [DJa) and using Eq. ([5]), we 
obtain Ag = 0.014. Unexpectedly in a doped sample 
with fast electron-hole recombination, the overall spin- 
relaxation rate increases with increasing pump power; 
see Fig. [DJc). In a high- mobility sample, a decrease in 
the spin relaxation-rate with increasing initial spin polar- 
ization has been observed, which goes into the opposite 
direction. '^^ In an attempt to minimize this pump-power 
dependence, we used a low pump power of 50 /iW for the 
measurement in Fig.ISjJa). 

Table [J summarizes the contributions to the 
anisotropic I/T2 for T = 30K. The calculated sum of the 
relaxation rate is by about 0.8 x 10^ s~^ lower than the 
measured value. This discrepancy indicates the presence 



of another spin-dephasing mechanism. A possible candi- 
date is a random SO field originating from the Coulomb 
potential of ionized dopants or from surface roughness 
of the QW. It has been pointed out that such spatial 
fluctuations might limit the spin lifetime for symmet- 
ric (110) QWs (Ref.lsil) or in the case where a = 
The importance of this effect is probably smaller in our 
samples where both a and /3 are finite but not equal 
in size. In a small-gap semiconductor, the Elliott- Yafet 
(EY) mechanism contributes to the spin dephasing»2£ By 
estimating the importance of the EY mcchanisroi^ in our 
sample using the measured mobility and the known band 
parameters, we obtain a spin-relaxation rate on the or- 
der of 5xl0^s^^. This is negligibly small, but there are 
indications'^^ that the EY spin-dephasing rate might be 
larger than estimated with the equation derived for a 
bulk semiconductor. 

The weak variation in I/T2* with T shown in Fig.[SI^a) 
can be understood as a consequence of little temperature 
dependence of the individual contributions to 1 /Tj* ■ As 
pointed out in Ref.Hil, the DP dephasing rate depends 
only weakly on T in the degenerate regime and in the 
intermediate temperature range, apart from its propor- 
tionality to the electron scattering time. As our mobility 
is quite constant over the temperature range measured, 
we do not expect large variations. We observe no evi- 
dence of a dependence of the g factor spread on T. For a 
degenerate electron density and a constant mobility, also 
the T dependence of the EY mechanism should be small. 
The observed weak T dependence is therefore not surpris- 
ing, and has also been observed in other experiments. ^'i 

To conclude, we have measured the SO interaction co- 
efficients a and /3 as a function of T and find no signif- 
icant T dependence. From a and /3, the measured Hall 
mobility and the anisotropy in I/T2, we estimate the 
contribution from DP spin dephasing, and find that DP 
alone cannot explain the measured I/T2*. From a linear 
increase in I/T2* with -Bextj we identify an inhomoge- 
neous broadening from a spread in the electron g fac- 
tor. These effects do not account for all of the measured 
spin-dephasing rate. We speculate that EY or an inho- 
mogeneous SO field might induce an additional isotropic 
contribution. A more detailed study of the nature of 
the elastic and inelastic electron-scattering mechanisms 
involved might facilitate an exact attribution to the dif- 
ferent decay mechanisms. 
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